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ABSTRACT 

We analyze simultaneously six composite RXTE spectra of GX 339-4 in the hard state comprising 77 
million counts collected over 196 ks. The source spectra are ordered by luminosity and span the range 
1.6% to 17% of the Eddington luminosity. Crucially, using our new tool pcacorr, we re-calibrate 
the data to a precision of 0.1%, an order of magnitude improvement over all earlier work. Using 
our advanced reflection model relxill, we target the strong features in the component of emission 
reflected from the disk, namely, the relativistically-broadened Fe K emission line, the Fe K edge and 
the Compton hump. We report results for two joint fits to the six spectra: For the first fit, we fix the 
spin parameter to its maximal value (a* = 0.998) and allow the inner disk radius R ln to vary. Results 
include (i) precise measurements of R m , with evidence that the disk becomes slightly truncated at a 
few percent of Eddington; and (ii) an order-of-magnitude swing with luminosity in the high energy 
cutoff, which reaches > 890 keV at our lowest luminosity. For the second fit, we make the standard 
assumption in estimating spin that the inner edge of the accretion disk is located at the innermost 
stable circular orbit (i? in = i?isco) and find a* = 0.951^05 ( 90% confidence, statistical). For both fits, 
and at the same level of statistical confidence, we estimate that the disk inclination is * = 48 ± 1 deg 
and that the Fe abundance is super-solar, Ap e = 5 ± 1. 


1. INTRODUCTION 

GX 339-4 is one of the most thoroughly studied of the 
roughly 50 known black-hole X-ray binaries. Its orbital 
period is around 1.7 days, and for the best candidate pe¬ 
riod of 1 .7557 days the mass function is 5.8 ± 0.5 M 0 
(Hynes et al. 2003). Like nearly all black hole binaries, 
the X-ray source is transient, having undergone more 
than a dozen o utburst cycles since it s discovery in the 
early 1970s by |Markert et al.| (|1973|). During a cycle, 
GX339-4 often exhibits all know n X-ray states, which 
unfold in the canonical patt ern (Remillard & McClin¬ 
tock 2006 Dunn et al. 2010). During the rising phase, 
the source can reach exceptional luminosities in the hard 
state, which is the focus of this paper. 

The hard state is strongly dominated by a hard power- 
law component (T ~ 1.6). The thermal component, 
which contributes < 20% of the 2-20 keV flux, is faint 
and cool ( kT < 0^ keV) com pared to the thermal state 
(Remillard & McClintock 2006). The Fe K line is a ubiq¬ 
uitous spectral feature. Strong variability is a hallmark of 
the hard state (rms power > 10% in the band 0.1-10 Hz), 
while QPOs may be either present or absent. The state 
is associated with the presence of an AU-scale steady jet, 
and clear correlations between th e rad io and X-ray inten¬ 
sities are observed ((Corbel et al.|2013). A major question 
for the hard state is the geometry of the corona: While 
there is significant evidence that the corona in the hard 
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state is compact, it is quite unclear whether, e.g., it is 
ADAF-like and quasi-sphe rical, hugs the disk or origi¬ 


nates at th e base of a jet (|Corbel et al.||2000| |Shidatsu 
et aLpOn i. 


1.1. Controversy over the Location of the Disk's Inner 

Edge 

In the thermal state, there is abundant evidence that 
the accretion disk is truncated near the innermost sta- 


ble circular orbit (ISCO) (e.g. 

Gierlinski & Done 

2004 

Steiner et al. 

2010 Penna et a 

• 

2010 Zhu et al. 

2012) 

The standard 

paradigm tor the 

taint hard state is that 


as the luminosity decreases the inner edge of the disk re¬ 
cedes from the ISCO, leaving a hot advection-dominated 


& Yi 1994 


accretion flow (ADAF) or other coronal flow (Narayan 
(Narayan & McClintock 2008). While there 
ai 


is good evidence that at very low lumi nosities the disk 
is grossly trun cated (for a review, see Narayan & Mc- 
Clintock[2008), the location of the inner edge relative to 


the ISCO tor luminosities in the range ~ 0.1 — 10% of 
Eddington is a hotly-debated topic. With GX 339-4 as 
a principal test bed, two methods have been widely used 
to estimate the radius R ln of the inner edge of the disk in 
the hard state: (1) modeling the component of emission 
reflected from the disk, principally the Fe K line; and 
(2) fitting the continuum spectrum of the accretion disk. 
The former method, which is addressed in the following 
section, is the central topic of this paper. 

Efforts to estimate the inner edge of the accretion disk 
in the low/harcl state via disk reflection go back far¬ 
ther, but the first strong indication that disks may re¬ 
main close to the ISCO in brig ht phase s of the low/hard 
state was made by Miller et al. (2006b). Based on fits to 
the thermal component, a number of papers claim that 
there is an optically thick disk that extends inward to 
the ISCO in the hard state ( Miller et al.|2006b|a Rykoff 
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et al. 2007 

Reis et al. 2009 

2010 

Reynolds et al. 2010). 

This claim is strongly contested 

by 

Done et al. (2007) 


questionable when o ne cons i ders that s elf-consistent disk 
coronal models (e.g., Steiner et al.|2009|) return larger val¬ 
ues of the inner-disk radius. More recently, Miller and 
coworkers have invoked extreme values of the spectral 
hardening facto r in making the case for an untruncated 
hard- state disk (|Reynolds & Miller 2013; |Salvesen et al. 
20131. This evidence lor the presence of such a disk does 
not appear to us compelling given the difficulties of ob¬ 
taining accurate estimates of R m by modeling a faint, 
cool (kT < 0.2 keV) thermal component that is strongly 
Comptonized and cut off by interstellar absorption. 

1.2. Reflection Spectroscopy 

The reflection spectrum results from the reprocessing 
of high-energy coronal photons in the optically-thick ac¬ 
cretion disk. The result is a rich spectrum of radiative 
recombination continua, absorption edges and fluores¬ 
cent lines, most notably the Fe K complex in the 6-8 
keV energy range. This reflected radiation leaves the 
disk carrying information on the physical composition 
and condition of the matter in the strong fields near the 
black hole. The Fe K emission line (and other fluorescent 
lines) are broadened and shaped by Doppler effects, light 
bending and gravitational redshift. By modeling the re¬ 
flection spectrum, one can estimate both the disk inclina¬ 
tion and the dimensionless spin parameter a* = cJ/GM 2 
(—1 < a* < 1). In measuring a*, one estimates the ra¬ 
dius of the inner edge of the accretion disk and identifies 
it with the radius of the innermost stable circular or- 
bit, Rtboo, which sim ply and monotonically maps to a* 
( Bardeen et al.|[l972 ). For the three canonical values of 
the spin parameter, a* = +1, 0 and -1, Risco = Id/, 
6M and 9 M (for c = G = 1), respectively. 

The reflection model most widely used in the past for 
both gener al application and me asuring black hole spin is 
reflionx ( Ross fe Fabian|2005 ). Recently, an improved 
reflection model has been developed, RELXILL R which is 

fa Sz Kallman 


based on the reflection code xillver (Garcia 

2010 Garcia et al.||2011 2013 2014a I, and the relativis¬ 

tic line-emission code RELLINE (Dauser et al.|2010} 2013' 


2014|. Compared to reflionx, RELXILL incorporates a 


superior treatment of radiative transfer and Compton re¬ 
distribution, and it allows for the angular dependence of 
the reflected spectrum. Furthermore, by implem entin g 
the routines of the ph otoionization code XSTAR (Kall¬ 
man & Bautista 2001), RELXILL provides an improved 
calculation of the ionization balance. At the same time, 
limitations of the model include assuming that the den¬ 
sity of the disk is independent of vertical height, that the 
illuminating radiation strikes the disk at a fixed angle of 
45 deg, and that apart from Fe all the elemental abun¬ 
dances are assumed to be solar. The results presented 
in this paper were derived using relxill to model the 
relativistically-blurred reflection component from the in¬ 
ner disk and xillver to model a distant reflector. 

It is important to appreciate the faintness of the re¬ 
flected features that are crucial for probing effects in the 
regime of strong gravity, the features that one relies on 
for estimating R ln and constraining black hole spin. For 

4 http://www.sternwarte.uni-erlangen.de/research/relxill 


example, in the spectrum of GX 339-4, even the most 
prominent feature, the Fe K line, has a typical equiva¬ 
lent width of ~ 0.1 keV, and the peak intensity of the 
line is only about 10% of the local continuum (Section[3]). 
Sensitivity to such faint features requires both high-count 
spectra and a well-calibrated detector. 

1.3. The Special Quality of This Study 

The principal detector aboard the Rossi X-ray Timing 
Explorer ( RXTE ) was the Proportional Counter Array 
(PGA), which was comprised of five nearly identical Pro¬ 
portional Counter Units (PCUs), each with an effective 
area of 1600 cm 2 and with sensitivity from 2-60 keV. 
Despite the limited spectral resolution of the instrument 
(ss 17% at 6 keV) the archive of PGA data amassed 
during the RXTE mission (1995-2012) continues to be 
preeminent for the synoptic study of stellar-mass black 
holes. A few-dozen bright black holes were observed daily 
during their outburst cycles with typical exposure times 
of a few ks. Some 15,000 individual spectra were ob¬ 
tained with a net total exposure time of 30 Ms (1 year). 
In this paper, we report the results of our analysis of six 
hard-state spectra of GX 339-4, each a summation of 
dozens of individual exposures (Section[3]). For the spec¬ 
trum obtained at maximum luminosity (i/^Edd = 17%) 
with an exposure time of 46 ks, the total number of 
counts is 40 million and the counts-per-keV in the con¬ 
tinuum at 6.4 keV is 4.4 million, while the total number 
of counts in the Fe K line region (3-10 keV) is 28 million. 

A limitation of the PCA, which has not allowed the im¬ 
plied statistical precision to be realized in modeling data, 
has been the ap precia b le ~ 1% uncertaint ies i n the de¬ 
tecto r response ( Jahoda et al.||2006 Shaposhnikov et al. 
2012|) . We have overcome this limitation by developing 
a calibration tool, called PCACORR, that increases the 
sensitivity of the RXTE PCA detector to faint spectral 
features - such as the Fe K line/edg e - by up to an order 
of magnitude (Garcia et al. 2014b). By applying PCA¬ 
CORR to a large number of spectra for three black holes, 
we found that the tool improved the quality of all the 
fits, and that the improvement was dramatic for spectra 
with > 10' counts. The tool allows one to achieve a pre¬ 
cision of ~ 0.1% rather than ~ 1%, thereby making full 
use of spectra of bright sources with ~ 10 6 counts per 
channel. 

Consequently, our study of the reflection spectrum of 
GX 339-4 greatly im proves o n ea rlier s tudies using the 
PCA, such as that by Plant et al. (2015). A limitation of 
PCA data is its modest resolution, while its major advan¬ 
tage is its freedom from the problematic effects of pileup, 
which is commonly a serious problem in analyzing and 
interpreting data for b right sources obtained using CCD 
detectors (see Section 6.1.3). Another advantage of the 
PCA, which has only recently been matched by NuSTAR, 
is its high-energy coverage, which allows observations of 
both the Fe K region and the Compton hump using a 
single detector. 

This paper is organized as follows: Section [2] describes 
the observations and data reduction, and Section [3] out¬ 
lines our procedure for combining the individual spectra 
into six composite spectra. The luminosities of these 
spectra, which we refer to throughout as Spectra A-F, 
range over an order of magnitude. Fitting the spectra in¬ 
dividually, while emphasizing the importance of correct- 
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X-ray Hardness 

Fig. 1. Hardness-intensity diagram for all PCU-2 RXTE observations of GX 339-4. The vertical axis shows the raw PCU-2 count rate 
(for reference, 1 Crab ~ 2600 counts s 1 ), a proxy for the X-ray intensity and luminosity. Plotted on the horizontal axis is th e hardness 
ratio HR defi ne d as the ratio of source counts at 8.6—18 keV to the counts at 5—8.6 keV. Following further the conventions of |Remillard| 
& McClintock (2006), the hard-state data considered exclusively in this paper are defined to have HR > 0.75. The six boxes labeled A-F 
define the data sets we sum to create Spectra A-F, which are used in our analysis throughout the paper. 


ing the data using the PCACORR tool, is the subject of 
Section [4 Our key results appear in Section [5] Therein, 
we describe how we fit Spectra A-F simultaneously, first 
fixing the spin parameter and letting the inner-disk ra¬ 
dius vary, and then allowing the spin parameter to vary 
while fixing the inner radius at the ISCO. We discuss our 
results in Section[6]and offer our conclusions in SectionO 


2. OBSERVATIONS AND DATA REDUCTION 

Our reduction and analysis of the R XTE PCA data 
are detailed in Ga rcia et al. (|2014b ), which follows 
the procedures of |McClintock et al. (|2006|). The 
data were obtained in “Standard 2” mode and seg¬ 
mented into contiguous intervals with exposure times 
ranging from 300 s to 5000 s. Background spec¬ 
tra, which were derived using pcabackest and the 
model pca_bkgd_cmvle_eMv20111129 .mdl, were sub¬ 
tracted from the data. Response files were generated 
using pcarmf (version 11.7) and the energy-t o-channel 
conversion table (v ersion e05v04) described in Shaposh- 
nikov et al. (2012). Throughout, we analyze just the 
data coliectedusing the best-calibrated detector, PCU- 
2, which also provides the richest data set. 

A s a cruci al final step, we apply the tool PCACORR 
(Garcia et al. 2014b) to the data and thereby calibrate 
the detector to a precision of ~ 0.1%; we include a sys¬ 
tematic error of this magnitude in all our data analysis. 
This ste p gre atly enhances the sensitivity of the detector 
(Section 4.31 to the reflection features that are our focus. 

Figure [T] presents a hardness-intensity diagram, PCA 
count rate versus PCA hardness ratio HR, for all 1471 
RXTE PCU-2 observations of GX 339-4 obtained be¬ 


tween 1996 July 26 and 2011 April 5. The hard-state 
data considered exclusively in this paper are defined to 
have HR > 0.75. In order to boost the signal-to-noise, 
we define the six boxes A-F shown in Figure[I] Each box 
contains a number of spectra, all of them corresponding 
to roughly the same source intensity. We combine all the 
spectra within a box using the procedures described in 
Section [3j importantly, we do not combine spectra ob¬ 
tained during different outburst cycles. Except for Box 
A, which is comprised of observations taken during the 
2002 outburst, all the other boxes contain observations 
taken during the 2010 outburst. Ultimately, we produce 
six master spectra (A F), one for each box. 


3. COMBINING SPECTRA 

We now outline our procedure for combining the indi¬ 
vidual spectra in a box to create Spectra A-F in such 
a way as to eliminate small variations in the power-law 
index under the assumption that the reflection features 
are unaffected by small changes in the continuum. For 
each box separately, we first fitted the individual spectra 
to a simple absorbed power-law (Tbabs*powerlaw) using 
a fixed hydrogen column density of IVh = 3 x 10 21 cm -2 , 
which is simil ar to the expected col umn in the direction 
of GX 339-4 (Kalberla et al. 20051. No evidence for a 
thermal component was found m any of the spectra. The 
fits were performed in the 3-45 keV band where the Fe 
K features are most pronouncec0 We then created in¬ 
dividual residual spectra (data counts minus model) and 
summed them, thereby greatly enhancing the residual 
features present in these spectra. 


5 Inclusion or exclusion of the 4—7 keV band was found to have 
a negligible effect on our final results. 
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Figure [2] shows the residual spectra for the six boxes. 
The striking features in each spectrum are the Fe K line 
and the K edge, which are revealed with precision in 
these high signal-to-noise spectra. Surprisingly, the over¬ 
all structure of the residuals are in all cases quite similar, 
despite the factor of ~ 10 spread in luminosity (Figure [I]). 
Upon closer examination, however, one sees that the line 
width, the position of the edge, and the shape of the 
Compton hump differ to some degree among th e bo xes. 
This point is discussed in more detail in Section [5T| 

We now use these residuals to create Spectra A-F, 
which constitute our prime data set. To be specific, con¬ 
sider the creation of Spectrum A: For box A we generate 
a single template continuum spectrum using the aver¬ 
age values of the fit parameters T and the normalization 
N. This spectrum is generated synthetically using the 
f akeit task in XSPEC; its net exposure time is the sum of 
the exposure times of all the 23 individual spectra in box 
A. Finally, we add this continuum spectrum to its cor¬ 
responding summed residual spectrum to complete the 
generation of Spectrum A. 

The resulting spectrum is superior to that obtained by 
simply summing the individual spectra directly because 
it seamlessly eliminates the effects of small differences in 
the power-law index and normalization among the spec¬ 
tra. The mean power-law parameters and other infor¬ 
mation describing these six spectra, which are hereafter 
our focus, are summarized in Table [T] We compute the 
luminosity for a spectrum using our model fluxes in the 
1-100 keV band and assuming a distance of D = 8 kpc 
and black hole mass of M = 10 M 0 (for details, see 
Footnote a to Table 1). 


4. FITTING SPECTRA A-F INDIVIDUALLY 

The residual plots in Figure [2] unambiguously demon¬ 
strate that a strong reflected component is present, which 
is widely attributed to the illumination of the disk by a 
hot corona. Invoking this paradigm, we proceed to fit 
each of the six spectra using our phys ically-m otivated 
reflection code relxill v0.2g (Garcia et al. 2014a). As 
before, Galactic absorption is modeled using Tbabs with 
fixed column density (ATh = 3 x 10 21 cm -2 ) . For the 
Tbabs model ( Wilms et al. 12000 ), we used the A nders fe 
Grevessel fl!9 89|) set of solar abundances and tnefVerner 
et al.| (|199()|) photoelectric cross sections. 

We fit Spectra A-F in turn to a succession of four mod¬ 
els; the final adopted model in each case yields a good fit 
with \ 2 /i/ ~ 1. Table fn provides detailed information on 
the quality of the fit for each spectrum and each model. 
In Figure [3j we show for Spectrum A with 4 x 10 7 counts 
the most challenging case - residual plots for the pro¬ 
gression of the four models, which we now describe. 

Model 0: Tbabs*powerlaw. An absorbed power-law 
model, which is clearly deficient, prominently displays 
the principal reflection features, the Fe K line/edge and 
Compton hump, in the residuals (Figure J3j). 

Model 1: Tbabs*relxill. A greatly improved fit 
to all six spectra is achieved by replacing the power-law 
with our fully relativistic reflection model. For simplicity 
and to achieve definiteness, we fix the spin to its extreme 
value of a* = 0.998 and assume the canonical dependence 
of disk emissivity with radius, namely oc r~ 3 . This model 
already delivers fits of reasonable quality (Table[2]) . Some 
pronounced residual features remain, which are most evi¬ 


dent for the most luminous case, Spectrum A (Figure pi). 
Specifically, two apparent absorption features flank the 
Fe K line at ~ 5.6 and ~ 7.2 keV. Other features are also 
present at higher energies in the region of the Compton 
hump (~ 20 — 45 keV). 

Model 2: Tbabs*(relxill+xillver). The residu¬ 
als are significantly reduced by including an unblurred 
reflection component via xillver. Physically, this re¬ 
flector could be cold material in a wind or in the outer 
region of a flared disk (see below). The xillver param¬ 
eters are linked to those of relxill with two exceptions: 
The ionization parameter was fixed at its minimum value, 
log £ = 0, and the Fe abundance was fixed to solar (i.e., 
Ap e = 1). Linking the Fe abundance results in a signif¬ 
icantly worse fit and a compromise value of abundance 
that is midway between the low value required by the 
unblurred component and the super-sol ar valu e required 
by the blurred component (see Section 6.1.4). This re¬ 
sult, and the uncertain origin of the unblurred compo¬ 
nent, motivate our choice of solar Fe abundance for the 
distant reflector. We have no good explanation for the 
different Fe abundances required in fitting the blurred 
and unblurred reflection components. Further discussion 
on the Fe abundance is presented in Section [6.1.4| 

Inclusion of the xillver component, which introduces 
only one new free parameter, namely its normalization, 
quite significantly improves the fit to all the spectra ex¬ 
cept Spectrum F, which has the fewest counts. While the 
xillver component improves the fit at low energies and 
in the region of the Compton hump .a strong residual 
feature remains at ~ 7.2 keV (Figure[3|. 

Model 3: Tbabs* (relxill+xillver) *gabs. We 
model the remaining residual feature near 7.2 keV phe¬ 
nomenologically as absorption using a single Gaussian. 
The addition of this component improves the fits sub¬ 
stantially for Spectra A-C, i.e., those with many counts, 
while it has only a marginal effect for Spectrum D and 
a negligible effect for Spectra E and F (Table [ 2 ]) . As ex¬ 
pected, its importance is greatest for Spectrum A where 
it completely eliminates the strong 7.2 keV residual fea¬ 
ture (Figure[3]) and produces a very good fit to this spec¬ 
trum, despite its extreme statistical precision (4 x 10 7 to¬ 
tal counts), with an allowance for systematic error of only 
0.1% (Section^. 

It is important to note that the inclusion of the 7.2 keV 
feature has a significant effect on some important model 
parameters. In particular, we find that including the 
Gaussian component (Model 3) increases the inclination 
and decreases the inner-disk radius by about 4 degrees 
and 20%, respectively, compared to excluding the com¬ 
ponent. The changes in the other fit parameters are rel¬ 
atively much smaller. We adopt Model 3 as our fidu¬ 
cial model for all six spectra, thereby assuming that the 
7.2 keV absorption feature has a physical origin. While 
it is beyond the scope of this work to establish a definite 
physical interpretation of the feature, we now briefly con¬ 
sider some plausible explanations. 


4.1. The 1.2 keV absorption feature and the efficacy of 
the xillver component 

We first consider the likely possibility that the absorp¬ 
tion feature is largely an artifact related to the uncertain 
energy resolution of the PC A. We then discuss the one 
plausible physical explanation for the feature known to 
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Fig. 2.— Residuals (data-minus-model) are computed for each box (defined in Figureml by subtracting an absorbed power-law fit to the 
individual spectra in the box. The extreme statistical precision results from summing millions of counts, ranging from ~ 28 million for Box 
A to ~ 3 million for Box F (3-10 keV). Note that the spectra are all scaled differently; e.g., the peak signal for the highest-luminosity Box 
A is ~ 15 times greater than for the lowest-luminosity Box F. Remarkably, the appearance of the spectra is quite insensitive to luminosity. 


TABLE 1 

Properties of Spectra A-F and the boxes defined in Figure [Tj 


Spectrum 

L/L e dd (%)“ 

Count Rate 

# Spectra 

< r > 

<N> 

Exp. (ks) 

A 

17.3 

1000 - 1100 

23 

1.72 ±0.01 

2.12 ±0.11 

45.7 

B 

14.2 

800 - 900 

7 

1.75 ±0.03 

1.81 ± 0.10 

10.3 

C 

11.9 

600 - 700 

11 

1.69 ±0.01 

1.22 ±0.06 

27.0 

D 

7.9 

350 - 400 

7 

1.61 ± 0.60 

0.60 ±0.03 

15.7 

E 

3.9 

150 - 200 

18 

1.52 ±0.01 

0.23 ±0.02 

24.6 

F 

1.6 

50 - 100 

43 

1.59 ± 0.08 

0.11 ± 0.03 

72.7 


“ Note. — Eddington-scaled luminosities assuming D = 8 kpc and M = 10 Mq (corresponding to I/Rdd = 1.25 X 10 39 erg s 1 ). and based 
on the fluxes computed over the 1 100 keV band using the model and fit parameters summarized in Table 3. 


us, namely that it is produced by absorption in a highly- 
ionized wind. Finally, we consider the role of the xillver 
component not only in modeling the residual features 
near the Fe K line, but also its role in improving the fit 
quite generally. 

4.1.1. On the accuracy of the PCA energy resolution 

The presence of residual absorption features bracket¬ 
ing the Fe line at ~ 5.6 keV and ~ 7.2 keV suggests 
the possibility that the PCA resolution may be better 
than assumed in generating the PCU-2 response. We 
have explored this possibility for Spectrum A. We test 
the effects of slight changes in the value assumed for the 
detector resolution by moderately smoothing the data, 
with the results shown in Figure [4j The smoothing is 
accomplished using a Gaussian kernel operating over the 
detector channels; the parameter / specifies the width of 
the Gaussian as a percent of the channel width. Accord¬ 
ingly, the curve in Figure [4] labeled / = 0 is unsmoothed, 
while the curves labeled / = 40 and / = 50 correspond 
to degrading the resolution of the data by 0.9% and 1.5% 
at 6.4 keV. 

Although this approximate approach to artificially tun¬ 
ing the detector resolution does not eliminate the residu¬ 


als flanking the Fe line, it does significantly reduce their 
strength. The test demonstrates that at this extreme 
level of statistical precision the fit to a line feature is 
very sensitive to the value assumed for the detector reso¬ 
lution. Specifically, if one assumes that the nominal value 
of resolution for the unsmoothed case (/ = 0) at 6 keV is 
17.0%, then the net resolution for / = 50 is 17.3% (i.e., 
the additional blurring has a width of one-half channel, 
equivalent to ~ 0.2 keV at 6.4 keV, which is combined in 
quadrature with the nominal resolution width). Mean¬ 
time, the resolution of the PCA is not known to sufficient 
accuracy to discriminate such fine differences (N. Sha- 
poshnikov, private communication). This suggests that 
the residuals near the Fe line may result from the PCA 
resolution being slightly better than assumed in modeling 
the detector response. However, this test is inconclusive. 
To properly assess the importance of tuning the resolu¬ 
tion, one must carry out a systematic analysis using the 
PCA calibration software, which is beyond the scope of 
this paper. 

4.1.2. On the possibility that the feature originates in a 
highly ionized wind 
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Fig. 3.— Data-to-model ratio (left) and contributions to x 2 (right) for Spectrum A resulting from fitting a sequence of four models. From 
top to bottom, the models increase in completeness and competence, starting with an absorbed power-law model to which, incrementally, 
is added a blurred reflection component (relxill); unblurred reflection (xillver); and a Gaussian absorption line gabs. The comparable 
residual plots for the other five spectra are qualitatively similar. 


TABLE 2 

Statistics of the individual fits to Spectra A-F 


Spectrum 

Model 


V 

xi 

Ax 2 /Au 

A 

0 

25094.90 

69 

363.694 



1 

299.07 

63 

4.747 

4132.64 


2 

151.07 

62 

2.437 

148.00 


3 

67.82 

60 

1.130 

41.63 

B 

0 

7653.27 

69 

110.917 



1 

117.00 

63 

1.857 

1256.05 


2 

81.16 

62 

1.309 

35.84 


3 

54.14 

60 

0.902 

13.51 

C 

0 

11849.29 

70 

169.276 



1 

142.44 

64 

2.226 

1951.14 


2 

93.21 

63 

1.480 

49.23 


3 

68.68 

61 

1.126 

12.27 

D 

0 

4880.44 

70 

69.721 



1 

89.91 

64 

1.405 

798.42 


2 

41.36 

63 

0.657 

48.55 


3 

35.69 

61 

0.585 

2.84 

E 

0 

2552.13 

70 

36.459 



1 

116.40 

65 

1.791 

487.15 


2 

65.97 

64 

1.021 

50.43 

F 

0 

2311.81 

70 

33.026 



1 

63.17 

65 

0.9719 

449.73 


2 

63.17 

64 

0.9870 

0.0 


If the 7.2 keV feature is not an instrumental artifact, a 
potentially plausible explanation is that it originates in 
a highly ionized wind that envelops the primary source. 
Disk winds have been observed in many black hole bi- 
nary systems, particularly at high accretion rates (e.g. , 


Ponti et al.|2012 Neilsen et al.|2012 Miller et al.|2006c|. 
We investigated this possibility by replacing the gabs 
component by the photoionized warm absorber model 


(warmabs). We forced the Fe Lya line at ~ 6.9 keV to 
be the dominant feature by setting the ionization pa¬ 
rameter to its maximum value (log £ = 5). We linked 
the Fe abundances of warmabs and relxill while the 
abundances of all the other elements remain at solar. 
The fitted blueshift of the Fe Lya required to model 
the 7.2 keV feature is z = 0.0576 ± 0.0101, which cor¬ 
responds to an outflow velocity of v = 1.7 x 10 4 km s -1 . 
The model provides a good fit (xt = 1-17), which is 
very comparable to that achieved using Model 3 (see 
the top-left panel of Figure [5] for details and a com¬ 
parison of the residuals). However, this interpretation 
seems unlikely on physical grounds due to the extreme 
column density required by the warm absorber, namely, 
iV|J bs = (7.7 ± 0.2) x 10 23 cm -2 . If one links the warmabs 
Fe abundance to that of the xillver component (i.e., 
Ap e = 1), the fit pegs at the hard limit of the warmabs 
model (10 24 cm -2 ). 

4.1.3. On the inclusion of the unblurred reflection 
component 

Our initial motivation for including the unblurred 
xillver component of reflection was the presence of the 
~ 5.6 keV and ~ 7.2 keV residual features flanking the 
Fe line. However, as Figure [3l makes clear, while the 
xillver component effectively eliminated the low-energy 
feature, it actually enhanced the 7.2 keV feature. The 
simplest ad hoc phenomenological approach to eliminat¬ 
ing both features is to include a pair of Gaussian ab¬ 
sorption lines in our model, which we did, fixing the 
widths of the Gaussians to 0.01 keV and allowing the 
energy and strength of each line to vary. While this 
model (Tbabs*relxill*gabs*gabs) does clean up the 
~ 5.6 keV and ~ 7.2 keV features, the quality of the fit, 
X 2 = 2.36 (x 2 = 139.25 for 59 d.o.f.), is much poorer 
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than that achieved with Model 3, xl = 1-13 (x 2 = 67.82 
for 60 d.o.f.), which uses xillver and a single Gaus¬ 
sian. As the top-right panel of Figure [5] makes clear for 
the stringent case of Spectrum A, thelatter model not 
only does a better job cleaning up the pair of targeted 
residual features, it is also more effective at improving 
the fit at most other energies as well. We consider this 
strong evidence for the presence in GX 339-4 of a distant 
reflector. 


4.2. Comparing combined versus summed spectra 

As discussed in Section [3j we combined the individ¬ 
ual spectra in order to increase our sensitivity to the 
reflection features while minimizing the effects of jitter 
in the power-law index. To verify the procedures we used 
in combining spectra, we compare the results of fitting 
our Spectrum A to those obtained by fitting a spectrum 
created by summing directly all the spectra in Box A[^] 
Applying PCACORR and fitting both spectra with our 
adopted model (i.e., Tbabs* (relxill+xillver) *gabs), 
we find that the model parameters are all consistent. 
However, the fit to the summed spectrum is of signifi¬ 
cantly lower quality (Ay 2 = 62.97) than the fit to Spec¬ 
trum A created using the procedures described in Sec¬ 
tion [3j Furthermore, as shown in the bottom-left panel 
of Figure [5] the residuals for the summed spectrum are 
larger in almost every energy channel. These results 
demonstrate that our method of combining the individ¬ 
ual spectra significantly improves the quality of the fit. 


4.3. The importance of applying the PCACORR tool 


As fully describ ed in Garcia et al. (2014b) and dis¬ 
cussed in Section L3 the PCACORR tool greatly re¬ 
duces the effects of instrumental features in PCA spec¬ 
tra, thereby making it possible to achieve good fits to 
high-count spectra at the 0.1% level of statistical preci¬ 
sion. The bottom-right panel in Figure [5] demonstrates 
the importance of applying the tool to Spectrum A, with 
its 4 x 10 7 counts. The figure compares residuals for a 
fit to uncorrected data to one using PCACORR-corrected 
data. To most clearly illustrate the power of PCACORR, 
we use Model 1 (Tbabs*relxill) and set the systematic 
errors to zero. Concerning the ~ 5.6 keV and ~ 7.2 keV 
features flanking the Fe K line, we note that they are 
present prior to the application of PCACORR, which con¬ 
firms that they are not introduced by the correction. The 
key message of the bottom-right panel in Figure [5] is the 
degree to which PCACORR diminishes these features and 
others, especially those below 10 keV and the one near 
30 keV, a feature that is likely related to the detector X e 
K-edge (Shaposhnikov et al. 2012 Garcia et al.||2014b ). 

5. FITTING SPECTRA A-F SIMULTANEOUSLY 


In order to achieve the strongest possible constraints 
on the key model parameters, we fitted Spectra A -F si¬ 
multaneously. This composite data set is an assemblage 
of 106 individual PCA/PCU-2 spectra of GX 339-4 in 
the hard state. The total number of counts is 77 million, 
34 million of which are in the 3-10 keV Fe K band. 

All fits are performed using Model 3 
(const*Tbabs*(relxill+xillver)*gabs). To ac¬ 
commodate the order-of-magnitude range of luminosity, 


6 The effects of changes in the detector response are negligible 
because the 23 observations were all made during a 10-day interval. 



Fig. 4.— Residual plots (data-model) illustrating the extreme 
sensitivity in fitting the Fe K line to the value assumed for the 
energy resolution of PCU-2. The residuals are for fits to Spectrum 
A using Model 1 (Tbabs*relxill). The resolution of the data have 
been slightly degraded by convolving them with a Gaussian whose 
smoothing width is characterized by the parameter f: The cases 
/ = 40 (green) and / = 50 (blue) correspond to decreases in the 
resolution of 0.9% and 1.5%, respectively, and / = 0 (red) is the 
unsmoothed case. 


we included a normalization constant that is unity for 
Spectrum A and floats for Spectra B-F. Where sensible, 
key physical parameters are tied: the spin a* of the black 
hole; inclination i of the system; line-of-sight column 
density Ah; Fe abundance Ap e ; and the normalization 
of the relativistic reflection component N r (see Table [3|. 
Given the uncertain origin of the absorption feature 
near 7.2 keV, which we model as a Gaussian, we also 
tie the central energy of this component while allowing 
its normalization to vary from spectrum to spectrum. 
The width of the Gaussian was fixed at 0.01 keV in all 
the spectra. Like the Gaussian normalization, all other 
model parameters are allowed to vary independently. 

Despite the extreme signal-to-noise of the composite 
spectrum, we must impose some additional assumptions 
in order to simultaneously constrain both the spin pa¬ 
rameter a* and the radius of the inner edge of the accre¬ 
tion disk i?i n . This is because these two quantities, which 
manifest almost indistinguishably in the red wing of the 
Fe K line profile, are extremely degenerate (Dauser et al. 


2013). Therefore, we have conducted two complement 


tary analyses, which we will refer to as JF-I and jf-ii, 
where JF signifies that these are joint fits (i.e., simulta¬ 
neous) to Spectra A-F (rather than the fits to Spectrum 
A alone that are featured in earlier sections). For JF-l, 
our aim is to constrain R- m , and we therefore keep the 
spin fixed at its maximum value of a = 0.998. For jf-ii, 
we tie the spin parameter for the six spectra and fit for 
it, while fixing the inner edge of the disk at the ISCO ; 
i.e., R\r, = Risen- We follow the guidelines of |Fabian| 


et al. (2012) and fix the emissivity index to its canonical 
value of 3; for both jf-1 and JF-II, we do not attempt 
to fit for this parameter because of the PCA’s limited 
spectral resolution. This choice is motivated by several 
JF-l and JF-II tests we performed that returned values of 
the emissivity index that were always < 4. 

The number of free fitting parameters is large, 52 for 
JF-l and 47 for JF-II. The complexity of the analysis 
dictated our approach: We performed Markov Chain 
Monte-Carlo (MCM C) runs using th e EMCEE-H AMMER 
Python package (|Foreman-Mackey et al. 20131, which 
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Fig. 5.— Comparative plot showing contributions to the total \ 2 (data-model) for each channel for fits to Spectrum A: ( top-left ) The red 
curve was computed for our adopted Model 3 (Tbabs* (relxill+xillver) *gabs) and is identical to the plot shown in the lower-right panel 
in Figure [3] The blue curve is for an alternative model that substitutes the warm absorber model wabs for gabs in our adopted Model 3. 
The alternative model provides a good fit to the data: x~ = 69.06 for 59 d.o.f. (xt = 1-130); compare y 2 values for Model 3 in Tables[3]and 
|4j ( top-right ) The red curve was computed for our adopted Model 3 and is identical to the plot shown in the lower-right panel in Figure [ 3 ] 
rhe blue curve is for an alternative model that replaces the unblurred (xillver) reflection component by a second Gaussian absorption line 
at ~ 5.6 keV (blue), which results in a distinctly inferior fit. As this residual plot shows, Model 3 performs better at almost all energies. 
( bottom-left ) Fits of our adopted Model 3 for two cases: (1) A direct sum of the 23 spectra in Box A (blue), and (2) a fit to Spectrum A 
(red), which was prepared by combining the spectra according to the procedures described in Section [ 3 ] While the model parameters are 
consistent for the two cases, Spectrum A provides a superior fit, as this comparison of the residuals makes clear. ( bottom-right ) Fits using 
Model 1 (Tbabs*relxill) for two cases: (1) The 23 individual spectra that comprise Spectrum A are corrected using the PCACORR tool 
(blue) and (2) they are left uncorrected (red). In this instance only, all systematic errors have been zeroed to most clearly illustrate the 
effect of applying PCACORR to the data. 


implements affine-invariant sampling. MCMC methods 
are powerful for high-dimensional analysis. Specifically, 
they enable an efficient exploration of parameter space 
and determine a posterior probability structure for the 
model of interest. 


5.1. Joint Fit I: Fixed spin and variable inner radius 

A principal goal of our study is to track the radius of 
the inner edge of the disk R m as the luminosity varies 
by an order of magnitude (i.e., over the range 1. 6-17 % 
of Eddington; Table [l])- As discussed in Section [lTT| a 
question of great interest is whether the inner disk is 
truncated in the hard state at low luminosities and, if 
so, to what extent. In order to be maximally sensitive 
to a disk that is only slightly truncated, we fix the spin 
to its maximum allowed value, namely a* = 0.998. In so 
doing, our focus is on determining how R ln trends with 
luminosity rather than obtaining accurate estimates of 
this parameter. We note that most spin determinations 
in the literature (which assume R ]u = i?isco) suggest 
that the spin is high (see Section 6.1.3). 

Figure [6] shows the fit residuals tor JF-I for two cases: 


(1) The top panel shows a data-to-model ratio for a fit to 
Model 0, i.e., the simple absorbed power-law model used 
to produce the residual plot for a fit to Spectrum A only 
(which is shown in the left-top panel of Figure [3]). This 
simple fit prominently displays the reflection features, 
which are strong for all six spectra. The profile of the 
Fe K line shows moderate variations among the spectra. 
As the luminosity increases, so does the intensity of the 
line (cf. Figure J2l). Additionally, the blue wing of the 
line extends to higher energies for the high luminosity 
spectra, which could be evidence for a shift in the Fe K 
edge caused by an increase in the ionization of the gas. 
There are obvious changes at high energies among the 
spectra, which are likely due to the evolution of the high- 
energy cutoff with luminosity. (2) Of chief interest, the 
middle and bottom panels of Figure [6] show the residuals 
for our adopted Model 3 (for a* = 0.998). 

Model 3 performs remarkably well for all spectra as 
indicated by the goodness of fit, xl == 1-06, and also by 
the uniformly ergodic appearance of the residuals across 
the energy band. The fit results are summarized in Ta¬ 
ble [3j A key result, which is discussed in detail in the 
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TABLE 3 

Results for jf-i: Fit parameters for Model 3, const*Tbabs*(relxill+xillver)*gabs, with fixed maximum spin a* and free Ri n . 


Model 

Parameter 

Spectrum A 

Spectrum B 

Spectrum C Spectrum D 

Spectrum E 

Spectrum F 

Tbabs 

relxill 

relxill 

relxill 

relxill 

gabs 

AT h (cm- 2 ) 
a* 

i (deg) 

^Fe 

N r 

E (keV) 



(7.0 ± 1.0) x 10 21 

0.998 

48.4 ±1.1 

4 6+ 0 - 5 

4 -°-0.3 

1 /1Q+0-05 

1-4o_o.o3 

7 l n+0.04 
'• iy -0.06 



Constant 

relxill 

relxill 

relxill 

relxill 

relxill 

xillver 

gabs 

Rin (7?ISCC>) 

r 

i°g£ 

E c ut 

R{ 

N x 

Strength 

i 

i 7 +O .2 

1.620 ±0.013 

q qi +0.03 
o.Ol—o.o7 

97±| 

0.21 ± 0.02 
n 97 + 0.02 

u - z ' -0.03 

0 021 +° 0° 9 
U.UZ 1 — Q 00 y 

0.91 ± 0.04 

1 5 +°' 4 
lo - 0.2 

1.682 ±0.016 

з. 24 ± 0.07 

129 ± 10 

0.22 ±0.03 

0.25 ±0.04 

0 029 +0 017 

и. uzy _ 0 012 

0.71 ±0.02 0.36 ±0.01 

i q+0.1 9 -1+0.3 

-L-O—0.2 z * A -0.5 

1.672 ±0.013 1.628 ±0.015 

o i 9+0.07 q r»qi +0.020 

°- iZ -0.03 o.UOi—o 013 

179 ± 14 660l^o 

0.21 ± 0.03 0.34tg;gg 

0.25 ± 0.04 0.34lg;gg 

0.036 ±0.016 0.05^0 03 

0.16 ± 0.01 
o 7+0.5 
+ '-1.5 

1.588 ± 0.010 

2.02 ±0.17 

> 840 

0.31 ± 0.03 

< 0.05 
°-° 8 to.05 

0.08 ± 0.01 
q 7 + 0.7 
‘ - 1.0 

1 C/lQ+0-007 

1.d4o_o 012 

2 05"*”°' 26 
z.uo_oio 

> 890 

0 'll ± 0-02 
u - ol -0.03 

< 0.03 

n 14± 0-10 
U ' i4 -0.05 


L/LEdd(%) 

17.3 

14.2 

11.9 7.9 

3.9 

1.6 


x 2 



402.49 




V 



379 




xl 



1.06 




Note. — For the given model components, the parameters from top to bottom are: hydrogen column density (-/Vh); dimensionless spin 
parameter (a* = cJ/GM 2 , where J is the angular momentum of the black hole); inclination of the inner disk (i); iron abundance with 
respect to its solar value (Ap e ); normalization of the blurred reflection component plus power-law continuum (AT r ); energy of the absorption 
Gaussian centroid ( E ); constant multiplicative factor between spectra; inner-disk radius (Rin), with Risco — 1.237 R g for a* = 0.998 
( R g = GM/c 2 ); power-law photon index (T); log of the ionization parameter (£ = 47 rF^/n, where F x is the ionizing flux and n is the gas 
density); high-energy cutoff (F C ut); reflection fraction (Rf, ratio of the reflected flux to that in the power-law, in the 20-40 keV band); 
normalization of the distant (unblurred) reflection (AT X ); strength of the absorption Gaussian; X-ray luminosity in terms of Eddington (see 
notes in Table [l| ; goodness of the fit (y 2 ); number of degrees of freedom goodness of the fit per degree of freedom (y 2 = X 2 / 1 ')- 
Uncertainties are based on a 90% confidence level. 


TABLE 4 

Results for jf-ii: Fit parameters for Model 3, const*Tbabs*(relxill+xillver)*gabs, with a* free and Ri n = Risco- 


Model 

Parameter 

Spectrum A 

Spectrum B 

Spectrum C 

Spectrum D 

Spectrum E 

Spectrum F 

Tbabs 

relxill 

relxill 

relxill 

relxill 

gabs 

7V h (cm -2 ) 
a* 

i (deg) 

^Fe 

Nr 

E (keV) 



(5.9t,;g) x 10 21 

0 95 + 0 ' 03 
u ' yo -0.05 

47.8+?'| 

5 4+ 1& 

°’ 4 -0.5 

1 44 + 0.04 

4-^—0.08 

7.23 ± 0.08 



Constant 

relxill 

relxill 

relxill 

relxill 

relxill 

xillver 

gabs 

Rin (7?ISC0) 

r 

log? 

Ecut 

Rf 

Wx 

Strength 

1 

1 604 +0 ' 010 
r.oud :_ 0 027 

3.33 ± 0.03 

92 tl 

0.20 ± 0.01 
r» 9q+0.05 
U-ZO — 0.02 

0 Q24+ 0 009 
U.UZ^4 _0 007 

0 90” t_u ' u ' 5 
u - yu -0.04 

1.658 ±0.018 

3.35 ±0.04 

118 ±8 

0 20”*” 0 02 
u - zu -o.oi 

0.26 ±0.04 

0.028t° ; °“ 

0 7i _l_u - u2 
J--0.03 

1 651 +0015 
1 - ool _ 0 . 0 22 

3-i6±g;Jg 

160 

0.20 ± 0.02 

0.24 ± 0.05 

0 Q37 + 0 016 
U.Uo / _q 012 

0.37 ±0.02 

1 

1 62”^ 002 
1 - oz -0.04 

3 05”*~® 04 
°- uo - 0.02 

440+ 230 
^ u -110 
n 97 +O.O 3 
u ‘ z/ -0.05 

0.13 ±0.06 

0 04+ 003 
d-u^_o 02 

0.17 ±0.01 

1 C7O+0.009 
1.0 < O—o 013 

1 96 +012 
1 - yo - 0.21 

> 830 

0 25”*" 0,05 
u,zo -0.03 

0.08 ± 0.06 

0.025 ±0.014 

0.08 ± 0.01 

1 cq7+0.009 
l.UO < _ 0>013 

2 0 +0 ' 2 
Z ' U - 0.2 

> 940 

n 90 + 0.02 

0 .^ 8 _o 04 

< 0.09 

0 020 +O ' 017 
u - uzu _ 0 . 009 


7'/iBdd(%) 

17.3 

14.2 

11.9 

7.9 

3.9 

1.6 


x 2 



418.66 




V 



384 




xl 



1.09 




Note. — For the given model components, the parameters from top to bottom are: hydrogen column density (./Vh); dimensionless spin 
parameter (a* = cJ/GM 2 , where J is the angular momentum of the black hole); inclination of the inner disk (i); iron abundance with 
respect to its solar value (Ap e ); normalization of the blurred reflection component plus power-law continuum (AT r ); energy of the absorption 
Gaussian centroid (F); constant multiplicative factor between spectra; inner-disk radius (Ri n ); power-law photon index (T); log of the 
ionization parameter (£ = AirF^/n, where F x is the ionizing flux and n is the gas density); high-energy cutoff (F cu t); reflection fraction 
(Rf, ratio of the reflected flux to that in the power-law, in the 20-40 keV band); normalization of the distant (unblurred) reflection (AT X ); 
strength of the absorption Gaussian; X-ray luminosity in terms of Eddington (see notes in Table [l|) ; goodness of the fit (y 2 ); number of 
degrees of freedom (y)\ goodness of the fit per degree of freedom (y 2 = x 2 /V). Uncertainties are based on a 90% confidence level. 
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Fig. 6 .— A simultaneous fit to Spectra A—F comprising a total 
of 106 individual spectra and 77 million counts, {top) Ratio plots 
for a joint fit (jf-i) using Model 0, the simple absorbed power-law 
model, which strikingly reveal the principal signatures of reflection. 
{middle) Ratio plots for these same data obtained by fitting our 
canonical Model 3 (Tbabs*(relxill+xillver) *gabs). This model 
produces an excellent fit with x 2 — 1-06. ( bottom ) Contributions 
to the total x 2 (data-model), again for Model 3 and the same data. 

following section, is evidence that the disk is moderately 
but significantly - truncated at the lowest luminosities 
that can be effectively explored using RXTE (i.e., at a 
few percent of Eddington). 


5.2. Joint Fit II: Constraining the spin of the black hole 


In order to obtain constraints on black hole spin using 
either leading method, reflection spectroscopy or contin¬ 
uum fitting, one must assume that Rj n = ^?isco (e.g., 


Reynolds 

2014 

McClintock et al. 

2014). In performing 

jf-ii, we make t 

his assumption for all six spectra in order 


to constrain the spin of the black hole. Doing so allows 
us to obtain a precise estimate of spin: a* = 0.95 +q 05 at 
90% confidence. As the summary of results in Table [4] 
shows, the other parameters are quite close to those ob¬ 
tained in JF-i (Table |3]), and the goodness of fit is of very 
comparable quality: y 2 = 1.09. Given the extreme sta¬ 
tistical precision, the residual spectra (data/model and 
y 2 ) for JF-II (not shown) are essentially indistinguish¬ 
able by eye from the spectra for JF-I, which are shown in 
Figure [ 6 } 

Concerning our spin estimate and the fundamental as¬ 
sumption that R'm = I?iscO; we again note that JF-I 
provides evidence for disk truncation at our lowest lu- 
minositietQ Meantime, these low-luminosity spectra are 
included in computing our single, tied JF -11 estimate of 
spin. The incorporation of disk truncation effects (which 
we have ignored) would imply an even higher spin value 
than is quoted above. 


6 . DISCUSSION 

The MCMC runs utilized 120 (jf-i) and 100 (JF- 11 ) 
“walkers,” each navigating a chain with a length of 

7 Furthermore, our JF-I results are consistent with some small 
degree of truncation even for our higher-luminosity data. 


100,000 elements, after having been initialized in a clus¬ 
ter distributed about the best fit. The first 50,000 el¬ 
ements of each walker were discarded in the “burn-in” 
phase during which the chain reaches its stationary state. 
The typical autocorrelation length, which is the interval 
over which the chain forgets its previous location, was 
several thousand elements; the corresponding net num¬ 
ber of independent samples of the parameter space was 
~ 10 4 . From the full distribution, we trivially obtain a 
probability distribution for any given set of parameters 
of interest by marginalizing over all the parameters that 
were outside that set. Flat priors were adopted for all 
model parameters. 


6.1. The four intrinsic parameters of the system 

We first discuss four parameters that are global for 
GX 339-4 (i.e., the same for all six spectra), namely, 
the Galactic hydrogen column density ./Vh, the spin pa¬ 
rameter a* of the black hole, the inclination i, and the 
Fe abundance Ap e . The entries for these (and all other) 
parameters estimated in the MCMC analysis given in Ta¬ 
bles^ and[4]are 90% minimum-width confidence intervals 
about the posterior maxima. 


6.1.1. Hydrogen column density 

Despite the limited low-energy coverage of the PGA, 
the hydrogen column density is well-constrained in both 


(Tables [ 5 ] 


jf-i and jf-ii to A h = (6.5_^ 5 ) x 10 : 
and [4]) , which is consistent with other estimates in the 
literature, including (in units of 10 21 


cm 


__ _ 4- 6 ( Kong 

et al. |2000l); 6 fl Zdziarski et al.|2004|); 5.4 (|Shidatsu et al. 
201 ip ; ana 5-8 ( Mendez fe Klis||1997 ). 


6.1.2. Inclination of the inner disk 

Previous estimates of inclination, which have been ob¬ 
tained by modeling the reflected component, run the 
gamut (see Table [5|) . Three papers using the same XMM- 
Newton and RXTE data uniformly report lo w values: 
Miller et al.| (|2006b|) and |Reis et al.| ([2008]) obtained 
i = 20 -?s deg and i < 20 deg, respectively, while |Done| 


& Diaz Trigo| (|2010|), using a different strategy for re- 
ducmg the data, found i ~ 20 — 27 deg. A much l arger 
inclination, i = 46 ± 8 degrees, was determined by |Shi- 
datsu et al. (2011) using Suzaku data. All of the se re¬ 


sults were obta ined using the reflionx mod els (Ross 
& Fabian 2005). Recently, Plant et al. (2015) fitted si¬ 


multaneously XMM-Newton and Suzaku data sets using 
both reflionx and xillver and reported two estimates 
of inclination: i = 36+@ deg and i = 42lg 4 deg. In a dif¬ 


ferent work, Plant et al. 


(2014b|) analyzed three low/hard 


state observations of GX 339-4 using a recent version of 
the relxill model and found i = 30+4 deg, a result 
that may be biased because it relies solely on low-energy 
XMM-Newton data. 

In JF-i, we obtained a tight constraint on inclination, 
i = 48.4 ± 1.1 deg (Table [3]), a result that is consistent 
with that obtained in JF-II (Table E]). Note that the 
relxill and xillver models used here properly treat 
the angular distribution of the reflected radiation, unlike 
earlier reflection models, w hich only provided an angle- 
averaged solution (Garcia et al)j |2014a I. 

As an aside, if one makes the usual assumption that 
the spin of the black hole is aligned with the orbital an- 
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46 47 48 

Inclination Angle 

Fig. 7.— A probability density map for the strength of the 
Gaussian absorption component (gabs) and the inclination angle 
from our MCMC analysis for jf-ii. While the inclination is tied 
between the six spectra A-F, the strength of the gabs component 
is free to vary. The distribution for each spectrum is coded by 
color. Only weak correlation is observed, with the gabs component 
becoming slightly stronger as the inclination increases and as the 
luminosity decreases. 


gular momentum vector (Fragos et al. 2010 Steiner & 
McClintock 2012), then the lower estimates of inclina¬ 
tion discussed above imply implausibly large values of 
black hole mass based on the Hynes et al. (2003) estimate 
of the mass function: e.g., M ~ 5.8/ sin^ i ~ 100 M 0 
for z ~ 25 deg. Meanwhile, our inclination implies 
M ru 15M 0 , consistent with the range of values observed 


for st ellar-mass black holes (Ozel et al. 12010 
20111 . 


Farr et al. 


The inclination angle is largely determined by the 
shape and position of the blue wing of the Fe K line, 
which in our fits is somewhat affected by the inclusion - 
or exclusion - of the Gaussian absorption feature. The 
energy of this feature was linked in all six spectra and 
constrained to 7.23 ± 0.08 keV (Table [3, while the nor¬ 
malization was free to vary. Figure 0 based on our 
MCMC analysis, shows that while the strength of the 
gabs component does increase by a factor ~ 2 with de¬ 
creasing luminosity, it only weakly interplays with incli¬ 
nation. This implies that this component, whose origin 
is uncertain (Section |4.1| ), has at most a modest affect 
on our estimate of inclination. 

6.1.3. Black hole spin 

The spin of the black hole has been estimated via 
reflection modeling using three independent data sets: 
a* = 0.939 ±0.004 using XMM-Newton/ EPIOMOS plus 
RXTE spe ctral hard-state data (Miller et al.|2006b |Reis 
et al. 12008 ); a* = 0.93±0.02 using XMM-X ewton/ EA’itd- 
pn plus RXTE sp ectral data in the very high (or steep 
power-law) state (Miller et al. 2004 Reis et al. 2008); 
and a* = 0.93 ± 0.01 (statistical) ± 0.04 (systematic) us- 
ing Suzaku data (Miller et al. 2008}. The corresponding 
estimates of inclination were all low (z ~ 10 — 20 deg) 
implying implausibly high estimates of black hole mass 
assuming spin-orbit alignment (Section 6.1.2). 

There is considerable uncertainty associated with these 
estimates of spin and inclination because of the effects of 
pileup, i.e., the arrival of two or more photons in the same 
or adjace nt CCD pixel within a sin gle frame time. For 
example, Done & Diaz Trigo (2010), analyzing precisely 


the same hard-state XMM-Newton/ EPIC-MOS data as 
Miller et al. (2006b) and Reis et al. (2008), conclude that 
the high spin reported by Miller et al. and Reis et al. is 
the result of severe pileup effects; using PN Tim in g-Mode 
dat a ( presum ably unaffected by pileup), Done & Diaz 
Trigo ( 20101 report evidence for a narrow Fe lin e and 
a truncated disk. Meantime, Miller et al. (2010) rebut 
the conclusions of Done and Diaz-Triga As a second 


example, the high spin reported by (Miller et al. 2008;) 
based on the ir analysis of Suzaku data (see above), is 
challenged by Yamada et al. (2009) who find - using the 
same data set - evidence tor a truncated disk and no 
n eed to invo ke a rapidly spinning black hole. 


Kolehmainen & Done (2010) applied the alternative 
continuum-fittingmethod to disk-dominated RXTE data 
collected during three different outbursts of GX 339- 
4. This method relies on accurate knowledge of the 
mass, distance, and inclination of the system, all of which 
are highly uncertain for GX 339-4. Using approximate 
bounds on these parameters, Kolehmainen & Done ob¬ 
tained “a strict upper limit” on the spin of a* < 0.9, 
which they claim is inconsistent with the spin estimates 
obtained by modeling the reflection spectrum. 

A chief virtue of the PGA data upon which we rely is 
its freedom from the confusing effects of pileup. A further 
virtue is the abundance of data, which allows us to track 
the behavior of GX 339-4 over a range of luminosity, as 
well as to reach extreme (~ 0.1%) levels of statistical 
precision. By assuming that the inner radiu s of the disk 
always remains at the ISCO (jf-II; Section 5.2), we es¬ 
tablished a firm constraint on the spin at a* = 0.951 qq5 
(90% confidence) while obtaining a precise estimate of 
the inclination, i = 47.8l°;4 degrees. Our spin result is 
in accord with the earlier Fe-line estimates, but our in¬ 
clination estimate is distin ctly d ifferent, and more in line 
with expectation (Section |6. 1.2). 

The results of our MCMC analysis allow us to search 
for possible degeneracies of the spin parameter with other 
fit parameters. Figure [8] shows for jf-ii spin probability 
distributions for three key parameters: (1) the inclina¬ 
tion angle, which affects the blue wing of the line; (2) the 
Fe abundance, which affects the strength of both the line 
and edge; and (3) the strength of the Gaussian absorp¬ 
tion, which also could affect the blue wing of the line. 
While there is no evidence for a substantial correlation 
between the spin and the inclination or the strength of 
the Gaussian, there is indication of a moderate positive 
correlation with the Fe abundance. From this positive 
correlation, it follows that an increase in Fe abundance 
will produce a diminished inner radius for f its p erformed 
with a fixed spin and variable i?i n (Section 6.3). 


6.1.4. Fe abundance 

The Fe abundance for the blurred reflection component 
(relxill) is surprisingly high: 4.6lg'3 and 5.dig'® in so¬ 
lar units for JF-I and JF-II, respectively (Tables [3| and |4|) . 
Most studies have merely assumed that the abundance 
is solar, while Allured et al. (2013) fitted Swift/RXTE 
hard-state data and afso found a super-solar abundance: 
Ap e = 2.4tJ;s2- The high abundance results directly 
from the remarkable strength of the Fe K line/edge rel¬ 
ative to the Compton hump (top panel of Figure [6]), 
which is usually the highest-amplitude feature in the re- 
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Fig. 8 .— Probability contours from our MCMC analysis of JF-II for the spin parameter and three other parameters of interest: inclination 
angle (left), Fe abundance (middle), and the strength of the Gaussian absorption component (gabs) for Spectrum A (right). A modest 
positive correlation is observed between the Fe abundance and the spin parameter. 


flection spectrum. A lower Fe abundance underpredicts 
the strength of the line/edge required to fit the Compton 
hump. Thus, our ability to constrain the Fe abundance 
is likely a consequence of our broad bandpass that pro¬ 
vides high-sensitivity coverage of all the principal reflec¬ 
tion features, from the Fe K line on through complete 
coverage of the Compton hump. This quality of cover¬ 
age is not provided by XMM-Newton data (even when 
RXTE data with a floating normalization are included) 
which may explain why others have not reported a super¬ 
solar Fe abundance. 

Quantitatively, forcing the Fe abundance to the so¬ 
lar value (Ap e = 1) results in a grossly unacceptable fit 
with xt = 9-9 and large residuals across the PCA band 
(Figure pi), while the inner-disk radius grows by about a 
factor of 10 compared to the fit with variable Fe abun¬ 
dance. We tried several alternative models (e.g., varying 
the emissivity index) in an unsuccessful attempt to find 
an acceptable model with lower Fe abundance. We note 
that fits to NuSTAR data for the most recent outburst of 


GX 339-4 likewise require a large Fe abundance (Fuerst 
eFak|M5 ). 



Energy (keV) 


Fig. 9.— Spectral fitting results, similar to JF-I, but here fixing 
the Fe abundance to the solar value (Ap e = 1). The top and bottom 
panels show respectively a ratio plot and contributions to \ 2 ■ This 
model fails to fit jointly the three principal reflection features (Fe 
K line and edge and the Compton hump), most noticeably for the 
higher-luminosity spectra, and the fit is unacceptable (xt = 9.9). 
A vertical dashed line marks the 6.4 keV rest-frame energy of the 
Fe K line. 


The Fe abundance has an influence on the shape of the 
Fe K line, and it may therefore in turn affect such param¬ 
eters of interest as the inclination or inner-disk radius. 
Moreover, the abundance also alters the continuum pho¬ 


toelectric opacity at higher energies, which modifies the 
depth of the Fe K edge and the red side of the Compton 
hump. These effects, which are subtle, are driving our 
fits because of the unprecedented signal-to-noise w e ha ve 
achieved. We shall return to this point in Section [6731 

We conclude that super-solar Fe abundance is a strong 
and inescapable requirement of these data. At the same 
time, the data also require the unblurred (distant) re¬ 
flection to have moderate - near-solar - Fe abundance. 
In our tests, fitting the spectrum with the most counts 
(Spectrum A), we found that adopting a single abun¬ 
dance for both blurred and unblurred components (i.e., 
relxill and xillver) results in Ap e ~ 3.6, an inter¬ 
mediate value between relxill (Ap e ~ 5), and the low 
value required for the xillver component. Though this 
approach may be intuitively more satisfying, it is strongly 
rejected by the data, with an increase in \ 2 of ~ 55. 
The effect of linking the abundances is to decrease the 
importance of the xillver component by reducing its 
normalization parameter from ~ 0.2 to ~ 0.07. An in¬ 
spection of the residual contributions to y 2 (lower panel 
of Figure [9]) reveals that the quality of the fit is degraded 
not only m the Fe K region but over the entire energy 
band. In summary, we find strong empirical evidence for 
the presence of an unblurred reflection component whose 
Fe abundance is much less than that of the blurred com¬ 
ponent. 

At the same time, there is no obvious reason why the 
inner-disk abundances should be so high. We note that 
similar physical processes may be occurring in AGN, 
since large Fe abundances are likewise found i n many 


cases when fitting relativistic reflection models fFabian 
2006L with 1H 0707-495 being a prime example (fltauseij 


et ai. 2012 Kara et al. 2015]) . Possibly, an unknown 
physical effect is being overlooked in current models that 
is artificially driving the Fe abunda nce to high values. 
For example, Reynolds et al. (2012) proposed that ra¬ 
diative levitation of Fe ions in the accretion disk at¬ 
mosphere could cause an apparent enhancement of their 
abundance. 


6.2. Parameters that evolve systematically with 
luminosity 

Setting aside the Gaussian absorption component (al¬ 
ready discussed in Section [4]) , for JF-I there are six im¬ 
portant parameters that are fitted separately for Spec¬ 
tra A-F: the inner-disk radius R ln ; the photon index T; 
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the ionization parameter £; the high-energy cutoff E cut ; 
the reflection fraction Rf, and the normalization of the 
unblurred reflection component (xillver). In this sec¬ 
tion, we show how these parameters depend on luminos¬ 
ity, and we discuss the causes of these dependencies. 

Figure 10 illustrates our MCMC results for JF-i (Ta- 
blef3J), the case of fixed spin. The probability distribution 
for each parameter is shown plotted versus the floating 
constant factor, which can be regarded as a proxy for the 
luminosity. The luminosity ranges over somewhat more 
than an order of magnitude. Each Spectrum is color- 
coded (see legend in top-left panel). The breadth of a 
distribution is a measure of uncertainty, while its shape 
indicates the degree of correlation of that particular pa¬ 
rameter with luminosity. We now discuss in turn the 
behavior of each parameter. 

6.3. Inner edge of the disk 

The evolution of the inner-disk radius R m with lumi¬ 
nosity is shown in the top-left panel of Figure [lOj Each 
spectrum delivers a good constraint on i?j n , allowing us 
to conclude that the inner edge of the disk moves out¬ 
ward by a factor of a few as the luminosity decreases by 
an order of magnitude, from a nominal value of 17% of 
Eddington to 1.6% of Eddington (Table [l]). 

This is a principal result of our paper because Ri n 
and its dependence on luminosity is a matter of central 
importance for the stu dy of black hole binaries in the 
hard state (Section 1.1). In Table [5] we summarize esti¬ 
mates of Rin in the literature for GX 339-4 in the hard 
state, while considering only those results obtained via 
reflection spectroscopy. The compilation includes results 
obtained using a wide variety of data and over a large 
range in luminosity (~ 0.1 — 20% of Eddington). At a 
glance, one notes the extreme range of values reported 
for Rin- The most notable conflict are two grossly dis¬ 
parate val ues rep orte d for the same XMM-Newton ob¬ 
servation: Reis et al. (2008) analyzed MOS and RXTE 
PCA data and reported i?j n = 2.041°'^ Ra. while |Plant 


be extremely truncated based on fits to a Suzaku and 
three XMM-NewtonRXTE spectra using xillver and 


et al. ( |2015[ ) analyzed EPIC-pn timing-mode data and 
reported R; n = 318±]| 5 R g . 

Figure |TT] shows all the values of Ri„ that appear in 
Table [5] plotted as a function of the Eddington-scaled lu¬ 
minosity. Several studies report results for multiple ob¬ 
servations over a range of luminosity; in these cases, the 
individual da ta points are highlighted in the left panel 
of Figure EH using colored tracks. Meanwhile, individ¬ 
ual measurements are shown in the right panel. Our 
results are shown in both panels with a solid red track 
connecting the data points. As noted above, we find 
that Ri n increases modestly with decreasing luminosity: 
Best-fit values trend upward from 2.1 R g to 4.6 R g as 
the luminosity decreases from 17% to 1.6% of Edding¬ 
ton. (Note that the values of Ri n in Table [3] are in units 
of Risco = 1.237R ff .) 

This trend is consistent with that found in previous 
studies except for that of Allured et al. (2013) (yel¬ 
low track), who fitted Swift and RXTE data using the 
reflionx model. Despite the general agreement that 
the inner radius shrinks with increasing luminosity, our 
estimates of R m at comparable values of luminosity are 
much smaller th an those repo r ted by others. 

For example, Plant et al. (2015) found the disk to 


et al. 

< 

ues 0 ] 



Kolehmainen 


radius (green track) by analyzing the same three 
XMM-Newton spectra (excluding the RXTE data) us- 
ing the rfxconv model (based on the reflionx tables; 


Kolehmainen et al. 2011). For their highest-luminosity 
data, they report two values of inner radius: One is a 
lower limit that is consistent with our results, R ln > 
2.2 R g , while the other (which includes a ~ 9 keV in¬ 
strumental feature in the fit; oran ge track ) is reason¬ 
ably consistent with the results of Plante^ al. (201 5]), 
Rin ~ 47 R„. Overall, the results of Kolehmainen et al. 


(20141 and Plant et al. (2015) arc similar, with the former 


authors reporting somewhat smaller values of Ri n . Inter¬ 
estingly, our results appear to be in reasonable agreement 
with an extrapo latio n of_the low-luminosity values of R; n 
reported by Petrucci et al. (2014]) (purple track), which 
is not necessarily expected since these authors assumed 
solar Fe abundance. 

Though the gross disparities in the reported values of 
R m may be partially due to differences in the models, 
this should be a secondary effect since, e.g., tests show 
that the models xillver and reflionx perform simi¬ 
larly (Garcia et al. 2013). The more likely reason for the 
inconsistent results is limitations of the data. One of the 
most severe of these is the effects of pileup, es pecially 
for the crucial XMM-Newton data (see Section 6.1.3). 
Another effect leading to major differences in results is 
whether or not high -ene rgy data were used. For exa mple 


the same EP 


Plant et al. (20151 and Kolehmainen et al. (2014) used 
G-pn data in timing mode, but while Plant | 


et al. (2015) used simultaneous RXTE data to extend the 


energy coverage, Kolehmainen et al. (2014) eschewed its 
use because of their concern over the cross-calibration of 
the two detectors. As a consequence of employing XMM- 


Newto n data only, the results of Kolehmainen et al. 


(20141 are highly sensitive to calibration issues associ¬ 
ated with the rapidly-falling and uncertain response of 
the EPIC-pn detector at energies > 9 keV. 

The PCA has important advantages despite its limited 
energy resolution and lack of coverage below 3 keV. Most 
notably, the PCA data are free from the contentious ef¬ 
fects of pileup that are inherent to CCD observations of 
bright sources. Meanwhile, the use of a single detector 
eliminates problems associated with cross-calibrating a 
pair of detectors. The much higher effective area of the 
PCA around the Fe line and Compton hump - and the 
many dozens of observations - yields spectra with orders- 
of-magnitude more counts than CCD spectra (Section^. 
Moreover, one can now fully utilize these many millions 
of counts per spectrum to detect subtle effects in reflec¬ 
tion features because the response of the PCA has been 
successfully calibrated to ~ 0.1% pre cision (Shaposh- 


nikov et al. 2012 Garcia et al. 2014b). These virtues 
of the PCA data are attested to by our success in fit¬ 
ting our reflection models (xt ~ 1) to six extremely 
high signal-to-noise spectra, which individually contain 
between 3 and 28 million total counts in the 3-10 keV 
band. Finally, the great abundance of data makes the 
PCA database unrivaled for synoptic studies of Galactic 
black holes. 

We now return to the question of the grossly discrepant 
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Fig. 10.— Variation of key model parameters with X-ray luminosity. The clouds of points in each panel (color-coded to correspond to a 
particular one of the six spectra) show the posterior density of the MCMC results for: the inner radius i?i n in units of the ISCO radius; the 
photon index T of the power law; the ionization parameter £; the high-energy cutoff E cu f, the reflection fraction Rf ; and the normalization 
N x of the unblurred reflection component xillver. The Constant Factor on the .r-axis. which is proportional to the Eddington-scaled 
luminosity, is normalized to unity (corresponding to L/L= 17 %) for Spectrum A (Table [TJ. 


results reported for R[ n (Table [5j Figure fTTl) while re¬ 
minding the reader that the Fe abundance affects the Fe 
K line profile and other reflection feature s at a detectable 
level given our signal to noise (Section 6.1.4). In turn, 
the Fe abundance affects other parameters, notably the 
inner-disk radius and spin parameter, which correlates 
positively with Fe abundance (Section 6.1.3). We now 
show that values of R m found by others are significantly 
biased by either the low signal-to-noise of their data or 
inadequate high-energy coverage. Such data make it dif¬ 
ficult to distinguish between small R[ n with large Ap e 
and la rge R ln with solar abundances, as we illustrate in 
Figure [12] which compares fits to two relxill models, 
one with Ap e = 1 and the other with Ap e = 5. The 
model with solar abundance (black curve) can only fit the 
data when the disk is strongly truncated, which serves 
to minimize the relativistic effects that blur the line pro¬ 
file. Note, however, that this model then fails to repro¬ 


duce the depth of the Fe K edge and underpredicts the 
continuum above ~ 30 keV. Figure [12] should be com¬ 
pared directly with Figure [9j where the limitations of 
the Ap e = 1 model are apparent from the fit residuals. 
Unlike most other data sets, the extreme signal in our 
data clearly discriminates between the two models. 

The truncation of the inner disk and the decrease in 
R m with increasing L/L^d, which we find, is a predic¬ 
tion of the advecti on-d ominated accretion flow (ADAF) 
model (see Section [O] ). In this paradigm, the inner disk 
evaporates becoming a very hot and optically thin ac- 


cretion flow that fills t he inner region (see, e.g., Meyer- 
Hofmeister et al.||2009). Our results are in line with this 


model, although our observations do not extend to the 
lower luminosities at which extreme truncation likely oc¬ 
curs. 

However, we note that these results are apparently at 
odds with our non detection of a thermal disk component 
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Fig. 11.— A comparison for GX 339-4 of our estimates with those in the literature (see Table |5|> of the inner-disk radius versus luminosity 
obtained by reflection modeling of hard-state spectra. The solid red track shown in both panels links our six measurements of Ri n (Table [M. 
(left) The variou s colored tracks sho w th e evolution of .Rin w ith Eddington-scaled luminosity (see footnote to Table m) reported by us and 
others (Tabic [ 5 ]). Plant ct al. (2015) and Kolehmainen et al. (2014) each present a pair of tracks, which are labeled CaJ and (b); the latter 
track in the case ofj Kolehmainen et al. (20141 includes a notch feature at ~ 9 keV. (rig ht) Single measurements of R\ n reported in a variety 
of studies (Table |5j. The track labeled Done+DiazlO corresponds to results reported by |Done & Diaz Trigo| ([2010| for fits to Epic-pn data. 
Data points with identical values of luminosity are slightly offset for clarity (precise values are given m TaDlejbJ. 



Fig. 12.— A comparison of two reflection models calculated 
using our relxill code: The red curve is for the parameters of jf-i 
(see Table [ 3 ]), in particular for Ap e = 5 and R[ n = Risco- The 
black curve is for a model with solar abundances (Ap e = 1) that 
has been doctored by increasing i^isco by 100-fold to best match 
the structure in the Fe K line. Meantime, the models are seen to 
differ greatly in the depth and structure of the Fe K edge and in 
the shape of the Compton hump. 

of emission. Making the usual assumption that all the 
observed power-law photons are generated by Compton 
up-scattering of disk photons, we would have expected 
to detect a thermal component for a disk that extends to 
such small radii (see Section [ 7 ]) . 


6.4. Parameters of the continuum: T and E cu t 

We find that the power-law photon index T is rela¬ 
tively constant despite the order of m agn itude increase 
in luminosity (top-right panel of Figure 10), a res ult that 
has been previously reported for GX 3$jF-4 (e.g. Wilms 


et al. 


1999; 
ige vai 


Zdziarski et al. 


ant et al. 20Ma). Its 


average value is 1.64U ± 0.035 tor jf-i and 1.625 ± 0.030 
for JF-II (std. dev., N = 6; Table [3| andjij) firmly i n the 
range for the har d state (1.4 < r < 2T; Remillard & 
McClintock||2006 ). 


In contrast to the constancy of the power-law index, 
the cutoff energy E cut systematically decreases with in¬ 
creasing luminosity from > 890 keV for Spectrum F down 


to 97 ±4 keV for Spectrum A (right-middle panel of Fig¬ 


ure 


10). This lower value of E cut for our highest lumi¬ 


nosity spectrum is of the same or der of magnitude as th e 


58.5 ± 2.2 keV value reported by Droulans et al. (2010), 
which is based on their analysis of simultaneous RX 7 E 
and INTEGRAL data obtained during another bright 
hard state of GX 339-4. 

Our model achieves good constraints for all six spectra. 
Remarkably, this is true for even the lowest-luminosity 
data (Spectrum F) for which the cutoff energy of > 890 
keV is far beyond the 45 keV limit of the PCA bandpass. 
This surprising result is a consequence of the detectable 
effects that are imprinted on the reflected component in 
the 3-45 keV band by photons with energies of hundreds 
of keV. We discuss the capability of the relxill model 
to probe the spectrum at extreme energies in Garcia et al.| 
(20151. 

in a Comptonized and isothermal corona, the high- 
energy cutoff is set by the electron temperature: E cut ~ 
(2 — 3 )kT e . In such a plasma, thermal disk photons 
are Compton up-scattered, thereby cooling the coronal 
electrons while producing the observed power-law con¬ 
tinuum. The slope of the power law depends on the in¬ 
terplay between the electron temperature and the optical 
depth r e , 


1 9 1 

2 + V 4 + d e T e ( 1 + T e / 3) ’ 


(1) 


(Lightman & Zdziarski,I1987), where d e = kT e /m e c 2 and 
m e c z = 511 keV is the electron rest mass. Values of 
these parameters for Spectra A- F, whic h are c onsistent 
with previous determinations (e.g. Wilms et al(p999]), are 
summarized in Table [6] for our nominal value"of the pho¬ 
ton index (T = 1.6). We find, as predicted by Equation]!] 
that the coronal temperature decreases with increasing 
luminosity, while the optical depth increases. 

6.5. The reflection fraction 

The reflection fraction Rf is a third parameter (in ad¬ 
dition to r and .Ecut) that provides information on the 
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TABLE 5 

Compilation of literature estimates of R in obtained by fitting reflection models to hard-state spectra of GX 339-4. 


Satellite Instrument L/L^d d (%) Ri n (Rg) i (deg) q High-energy? Ref. 


XMM-Newton EPIC-pn (TM) a 

1.42 

684t^ 

421/ 

3 

Yes 

[1] 



9721^3 

3611 

3 

Yes 

[2] 



1 in+ 80 
11U -40 

60 

3 

No 

[3] 



150*_ 50 

60 

3 

No 

[4] 

EPIC-MOS 

3.25 

5 ±0.5 

2D - *" 5 

zu -10 

3 

Yes 

[5] 



b- (N 
O O 

6 o 
+ 1 

o 

CN 

20_i.3 

3.16 ±0.5 

Yes 

[6] 

EPIC-pn 


10 ±2 

27 ± 3 

3 

Yes 

[7] 



60” 1 ” 40 

DU -20 

60 

3 

Yes 

[8] 

EPIC-pn (TM) a 


3181)® 5 

421/ 

3 

Yes 

[1] 



125±H 

3611 

3 

Yes 

[2] 



8911 

60 

3 

No 

[3] 



12811 

60 

3 

No 

[4] 


10.2 

15511 9 

421/ 

3 

Yes 

[1] 



79+42 

<Z -2l 

3611 

3 

Yes 

[2] 



2.2* 

60 

3 

No 

[3] 



47l)° 

60 

3 

No 

[4] 


< 0.05 

2111 

soil 

3 

No 

[14] 


< 0.05 

27l® 

3oil 

3 

No 

[14] 


< 0.05 

1611 

30l| 

3 

No 

[14] 

Suzaku XIS0,1,3/PIN 

0.14 

> 65 

18 

2-3 

Yes 

[9] 



> 798 

421/ 

3 

Yes 

[1] 



> 745 

3611 

3 

Yes 

[2] 



190 +l7u 

ib,u -90 

50 

2.3 

Yes 

[10] 



> 180 

20 

3 

Yes 

[11] 


0.13 

> 30 

20 

3 

Yes 

[11] 


0.19 

> 10 

20 

3 

Yes 

[11] 


0.25 

> 70 

20 

3 

Yes 

[11] 


0.60 

09 + 33 
oz -19 

20 

3 

Yes 

[11] 


0.91 

' * u -1.3 

20 

3 

Yes 

[11] 


2.0 

13 3 +6,4 
±o ‘°-6.0 

46 ± 8 

2.3 ±0.1 

Yes 

[10] 

Swift XRT 

1.33 

3-eli:S 

20 

3.2 ± 0.6 

Yes 

[12] 



6.7±l\ 

20 

3 

Yes 

[13] 


0.46 

< 10 

20 

3.1 ± 0.4 

Yes 

[12] 



3-6lS;g 

20 

3 

Yes 

[13] 


1.08 

12.8ll 0 8 

20 

3 

Yes 

[13] 


1.12 

6-9111 

20 

3 

Yes 

[13] 


1.68 

19.5ll 5 

20 

3 

Yes 

[13] 


2.23 

16-311/ 

20 

3 

Yes 

[13] 


5.21 

19.71/s 1 

20 

3 

Yes 

[13] 


“ TM=Timing Mode 

Note. — flllPlant et ah ( |2015| l implementing xillver; [2] IPlant et al. j _ __ 

[41 Kolchmaincn ct al. (2014) including a notch feature at 9 kcV; 5| Miller et al. (2006b); [61 Reis et al. 


imp lemen ting reflionx; 


31 Kolehmainen et al. (201 4|. 


201 U b : 8 [T )oik & Diaz Trigo ( 2 010b with fixed incl inati on; [91 [T omsick ct al. (2009); [10 Shidatsu et a I. I 120 ITT 
12 ] I'lomsick et ak] (|2UU»||; [13] Allured et al.| (|2013|) ; [14]|Plant et al.| ||2U14b|l. 


20081); l/l | Done &s Diaz Itig ol 
i~[l 1] jPetrucci et al7(i'20I4); 


TABLE 6 

Coronal properties 11 


Box 

£/T-Edd 
(%) 

E cut 

(keV) 

T 

(10 9 K) 

T e 

A 

17.3 

97 

0.45 

3.03 

B 

14.2 

129 

0.60 

2.50 

C 

11.9 

179 

0.83 

1.99 

D 

7.9 

660 

3.06 

0.72 

E 

3.9 

840 

3.90 

0.59 

F 

1.6 

890 

4.13 

0.56 


“ Note. - Assumes kT = l Lui and P = 1.6 (See Equation |l|). 


structure of the corona. In relxill, the parameter is 
empirically defined as the ratio of the reflected flux to 
the power-law flux in the 20-40 keV band. The re¬ 
sults of JF-I show that the reflection fraction of the rel- 
ativistically blurred component (relxill) ranges from 
0.2 < Rf < 0.3, decreasing modestly with increasing lu¬ 
minosity (bottom-left panel, Figure 10). This trend is 
surprising given that, at the same time, I?i n is decreas¬ 
ing so that the area of the reflector should increase. As 
a further w rinkle, one e xpects R{ > 1 based on simple 
arguments (|Dauser et al.(2014). 

There are several scenarios that can plausibly account 
for values of I?f < 1. We mention four and then discuss 
a new, alternative explanation. (1) An obvious expla- 
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nation is a severely truncated disk (more specifically, a 
disk with truncation radius large compared to the size 
of the corona). We discard this possibility as inconsis¬ 
tent with the small values we find for R ln (Section 6.3). 
(2) Another option is for the corona to be continuously 
outflowing at relativistic speeds, bea ming the bul k of its 
emission away from the disk (e.g. Miller et al. 2014} [Keck 
et al. [2015). While this is a possible explanation for 
R{ < 1, one outcome of this scenario is a resultant low 
value of the emissivity index, which is not obviously re¬ 
quired by our data. (3) The value of Rf may be depressed 
by our assumption of a cons tant-density disk atmosphere, 
as |Ballantyne et al. (2001) have shown in their studies 


of hydrostatic atmospheres. The hotter gas layer at the 
surface of a hydrostatic atmosphere additio nally scatters 


and blurs reflection features (see also Nayakshin & Kali- 
man 2001) thereby diluting the reflection signal relative 


to a constant density model. (4) The apparent strength 
of reflection features may also be reduced by the Comp- 
tonization of these features in an extended corona, as 


Wilkins & Gallo (20141 recently proposed. However, for 
such a corona to be effective in reducing i?f appreciably, 
it must have a large covering fraction which may inter¬ 
fere with detection of blurred reflection features from the 
inner disk. 

We propose an alternative explanation for Rf < 1 
based on the strong dependence of the reflected spectrum 
on the angle at which an illuminating photon strikes the 
disk. This angle crucially determines the characteristic 
depth in the disk at which the photon interacts; this in 
turn affects the limb - darkening/brighteni ng of the disk 
( Svoboda et al.|[2009| Garcia et al.||2014a ). A deficiency 
of reflionx, relxill, xillver and other widely-used 
reflection models is the simplifying assumption of a fixed 
incidence angle of 45 deg. However, a larger angle of in¬ 
cidence (measured with respect to the normal to the disk 
plane), for example, results in a hotter surface layer and 
therefore a weaker reflection signature (see Figure 5 in 
Dauser et al.||2013 ). 

'I’o test whether the assumption of near-grazing illu¬ 
mination substantially increases fitted values of Rf, we 
produced a new table of xillver reflection models with 
a fixed incidence angle of 85 deg and merged them with 
RELLINE to create a n ew high-incidence-angle version of 
relxill (Section 1.2). Fitting Spectra A-F as in Sec- 
(i.e., JF-l), the fit is slightly worse (Ay 2 = 9.71) 
istically comparable and still quite reasonable 
1.09). Notably, for the 85-deg model we find 
that Rf increases with luminosity and that i?f > 1 for 
the three most luminous spectra (A-C). Meanwhile, all 
the other parameters are consistent with those for jf-I 
(Table p} . Importantly, the Fe abundance remains un¬ 
changed!^ 

Figure [l3| compares the reflection factors computed for 
the two models. The large-angle model is more in accord 
with expectation, namely, the reflection fraction trends 
upward with luminosity and the values at the higher 
luminosities (with Rf n near the ISCO and with corre¬ 
spondingly large reflector area) are sensibly > 1. Thus, 
our results qualitatively suggest that the accretion disk 
in GX 339-4 is illuminated at near-grazing angles with 
respect to the surface of the disk. 

Within a few gravitational radii of the horizon, the 


tion 5.1 
but sta' 

(xl = 


1.6 

1.4 

1.2 

1.0 

0.8 

0.6 

0.4 

0.2 

0.0 


45 deg Incidence 
85 deg Incidence 


. I- 


S . 5- 


0.04 0.08 0.12 0.16 

Luminosity (LEdd) 


0.2 


Fig. 13.— Comparison of the reflection fraction versus lumi¬ 
nosity for two fixed value of the incident angle of the illuminating 
radiation: the 45 deg value (blue points), which is widely assumed 
in reflection modeling, and 85 deg (red points), the value we have 
assumed for this test. 

extreme bending of light rays causes pho tons to strike the 
disk over a wide range of angles (D auser et al.||2013 see 
in particular the middle panel of their Figure 5). Given 
the strong dependence of Rf on the angle of incidence, 
reaching firm conclusions concerning the reflection factor 
will require building a new generation of models, a task 
beyond the scope of this paper that will be addressed in 
future work. 

6 .6. Ionization parameter and geometry 

As expected, both the normalization N r and the ion¬ 
ization parameter £ of the blurred reflection compo nent 
(relxill) increase with luminosity (Table [3j Figure 10 
middle-left panel). In particular, the ionization param¬ 
eter changes from £ = 112.2 to £ = 2041.7, which 
traces very well the ten-fold increase in luminosity, from 
1.6% AEdd to 17% LEdd- For the strongly-illuminated 
portion of the disk, the variations in £ and L deviate 
mildly from the simple relation £ = L/nD 2 , where n is 
the density of the gas (fixed at n = 10 15 cm -3 for the 
relxill and xillver models used here[^J and D is the 
distance from the coronal source to the strongly heated 
portion of the disk. Thus, D increases only modestly 
as the luminosity decreases by an order of magnitude in 
passing from Spectrum A to Spectrum F: 


Df 

~D~a 


' L f ^a 
La£,f 


1.3. 


( 2 ) 


This small change is reasonable given the correspond¬ 
ingly mild increase in the inner radius obtained for JF-I: 


Rf 

Ra 


= 2 . 2 . 


(3) 


Likewise, the normalization N x of the unblurred re¬ 
flection com ponent (xillver) increases with luminosity 
(Figure [lOj bottom-right panel); although presumably 
the ionization parameter of this component also increases 
with luminosity, we approximate the state of the gas in 
the distant reflector as cold and neutral. 

8 The choice of gas density is relatively unimportant; it is the 
ionization parameter t hat largely determines the properties of th e 
reflected spectrum (see|Garcla &; Kallman|2010||Garcla et al.]2013|). 
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7. SUMMARY AND CONCLUSIONS 

We have presented an analysis of six composite RXTE 
PCA spectra of the X-ray binary black hole GX 339-4. 
All these spectra were taken when the source was in the 
hard state. The spectra correspond to luminosities rang¬ 
ing from 1.6% to 17% of the Eddington luminosity. The 
six spectra, each spanning the energy range 3-45 keV, 
comprise in total 77 million counts and a total exposure 
time of 196 ks. A unique feature of this work is our use 
of the tool PCACORR, which allows us to calibrate the 
PCA data to a precision of 0.1%. 

The spectra individually, and jointly, are well fitted 
by a model with three principal components: relxill, 
our model of relativistic ionized reflection; xillver, a 
minor component that models the effects of a cold, dis¬ 
tant reflector; and Tbabs, a standard model of Galactic 
absorption. We include an ad hoc Gaussian component 
(gabs) to model an absorption feature near 7 keV. The 
origin of this feature is unclear, but it is likely an ar¬ 
tifact resulting from a misestimate of the PCA energy 
resolution. 

We performed two joint fits of the six spectra. In the 
first of these, we fixed the spin to its maximal value, 
which allows the inner-disk radius R\ n to approach the 
IS CO radius, and we derived precise estimates for the 
evolution of R ln with luminosity. We find that the 
disk becomes increasingly truncated with decreasing lu¬ 
minosity. Specifically, as the luminosity ranges from 
17% to 1.6% of Eddington, R m increases from 2.1 R g to 


,(e.g., 

Petrucci et al. 2014 

Kolehmainen et al. 12014 

Plant 

et al. 

20151), our values ol K[ n lor comparable values ol 


The grossest disc repancy is the hundredfold larger values 


reported by Plant et al. (2015). 

That we rind such small values of the inner-disk ra¬ 
dius and no evidence for a thermal disk component is at 
odds with the current models. This is particularly true 
for Spectrum A with R ln = 2.1 R g . One expects such 
a modestly-truncated disk to be sufficiently hot (partic¬ 
ularly because it is heated by the corona; e.g., |Haardt| 
& Maraschi |1993 ) that we should have detected it with 
the PGA. This implies that our model somewhat under¬ 
estimates the true value of i?i n and that our model is 
incomplete. To address this problem, we are in the pro¬ 
cess of exploring an extended model that self-consistently 
treats the thermal, power-law and reflected components. 
This is a challenging problem whose solution is beyond 
the scope of this paper. 

Our analysis indicates that the factor of ~ 100 range 
in the values of R- ]n at fixed luminosity, which have been 
reported in the literature, is unlikely to result from the 
use of different reflection models; the shifts in R ln at¬ 
tributable to this cause appear to be relatively minor. 
Instead, the large disparity appears to be attributable to 
limitations of the data, one of which is the well-known ef¬ 
fects of pileup. In this paper, we highlight a particularly 
important effect, namely, the modest statistical quality 
of most data, which has resulted in observers fitting the 
blurred reflection component assuming that the Fe abun¬ 
dance is solar, whereas we demonstrate that super-solar 
Fe abundance is required for fits to data with extreme 
statistical precision. Specifically, we strongly constrain 


the Fe abundance (in solar units) to be Ap e = 5.0 ^q 4, 
which is the average value for our two joint fits. This 
strict requirement of the data is a promising and likely 
explanation for why, at luminosities ~ 1% of Edding¬ 
ton, we find evidence for relatively mild disk truncation 
compared to earlier studies. 

We acknowledge that the accuracy of our results are 
limited, systematically, by the pre senc e of an absorp¬ 
tion feature near 7.2 keV (Section 4.1) whose origin is 
unknown. However, our principal conclusions regarding 
the inner radius of the disk are sound, being subject to 
a minor uncertainty of about 20% arising from whether 
or not this feature is included in the model. 

As the source luminosity and the radiation field 
bathing the disk grow, the disk becomes increasingly 
ionized and its structure changes as R ln shrinks. At 
the same time, the large and steady decrease in the 
high-energy cutoff indicates that the illuminating coro¬ 
nal source is likewise evolving, as its temperature drops 
and its optical depth increases. 

In the second of our two joint fits to the six spectra, we 
made the standard assumption used in estimating black 
hole spin, namely, we fixed R ln to the radius of the ISCO. 
Doing so, we constrained the spin of the black hole to be 
a* = 0.95^o 05. We were able to achieve this statistical 
precision despite the limited spectral resolution of the 
PCA because of the quality of the data and its calibra¬ 
tion. 

If there is some truncation of the inner disk (i.e., 
i?i n > -Risco), then the spin is greater than the esti¬ 
mate given above. Our estimate of spin agrees well with 


(Miller et al. 2004 

2006b 

2008 

Reis et al. 

2008 

). It is, 

however, inconsistent with the upper limit o 

a * < 0.9 de- 

rived using the continuum-fitting method ( 

Oolehmainen 


curate values of black hole mass, disk inclination and 
distance that are required for successfully applying the 
continuum-fitting method are unknown for GX 339-4. 
Our result is also formally incomp atible with t he value 
of spin predicted for GX 339-4 by |Steiner et al.| (|2013 1 
based on the rel ationshi p between s pin and jet power 
proposed by Naray an fc McClintockl ( 2012), which has 
been challenged by|Russell et aT7| (|2()13|T 
We also obtain a precise estimate for the inclination of 
the inner disk of i = 48.I+4 3 deg. This value is further 
subject to an estimated systematic uncertainty of about 
4 deg arising from whether or not one chooses to include 
the 7.2 keV absorption feature in the model. Our value 
is inconsistent with the low values found earlier using 


the Fe-line method (Miller et al. 2004 2006b 2008 Reis 


et al. 

2008 

), while it is more 111 line with reasonable ex- 

pecta 

tions 

for the mass of the black hole based on the 


value of the mass function. 
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